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We formulate a theory of slow polaritons in atomic gases 
and apply it to the slowing down, storing, and redirecting of 
laser pulses in an EIT medium. The normal modes of the cou- 
pled matter and radiation are determined through a full di- 
agonalization of the dissipationless Hamiltonian. Away from 
the EIT resonance where the polaritons acquire an excited- 
state contribution, lifetimes are introduced as a secondary 
step. With detuning included various four-wave mixing pos- 
sibilities are analyzed. We investigate specifically the possi- 
bility of reverting a stopped polariton by reversing the control 
beam. 

PACS numbers: 42.50.Gy, 32.70.Jz, 42.50.Fx, 03.75.Fi 



Recently, electromagnetically induced transparency 
(EIT) 0-0] was shown to slow down dramatically M , or 
even to stop completely tJL[7|, laser pulses in atomic gases. 
The experiments involve media of three-level atoms in- 
teracting with two lasers — a control beam and a probe 
pulse. The atoms have two hyperfine ground states, \g) 
and | q), and an electronically excited state |e), as illus- 
trated in Fig. 1(a). Level g is populated initially, before 
applying the probe pulse to couple g and e. The role of 
the control beam is to introduce a transparency window 
so that the probe pulse propagates slowly in the medium. 

Such behavior can be understood in terms of a branch 
of slow polaritons appearing between two close atomic 
resonances, as considered by Juzeliunas (see Fig. 2(b) 
in H). Indeed, the control laser couples the states \q) 
and |e) dynamically, bringing level q into resonance with 
the excited level e. The excited level splits, then, into 
the doublet shown in Fig. 1(b), giving precisely the level 
structure required to form a branch of slow polaritons. 
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FIG. 1. (a) Atomic level scheme in the A configuration of 
EIT. (b) Equivalent scheme incorporating the dynamical cou- 
pling of q and e characterized by the Rabi frequency fi c . 



state physics. Over the last decade the polariton idea has 
been applied widely to describe the quantized radiation 
field in dielectric media [p-14|. Most studies, however, 
considered media of two-level atoms, and hence cannot 
accommodate the slow EIT polaritons. Slow polaritons 
appear in the analysis beyond two levels [p[ p^Jl4[ |; al- 
though the existing theoretical work does not deal with 
EIT, specifically. EIT (dark state) polaritons were first 
considered theoretically by Mazets and Matisov [g5| , and 
later by Fleischhauer and Lukin (ljj,[l7|] who suggested 
storing the probe pulse (stopping the polariton) by adi- 
abatically switching off the control laser. 

In this paper we present a systematic description of 
slow polaritons in EIT media. The theory is developed 
for atomic Bose-Einstein condensates (BECs) , but is also 
applicable to ordinary atomic gases. In contrast to previ- 
ous work |R5|-[r7|| , we explicitly diagonalize the full Hamil- 
tonian, including detuning from the EIT resonance and 
the contact interaction, and without making the rotating 
wave approximation in the interaction with the probe 
field. Away from the EIT resonance the polaritons ac- 
quire an excited-state contribution which leads to a finite 
radiative lifetime. With detuning included various four- 
wave mixing possibilities are suggested and we specifi- 
cally propose a scheme for reverting a stopped polariton. 
We also apply our formulation to the recent experiments 
MM where it is apparent from the explicitly constructed 
polariton modes that the storage and retrieval of the 
probe pulse takes place at the medium boundaries (to 
within an extremely small correction) and stopping the 
polariton can be achieved even by a sudden "turn-off" of 
the control laser. 

Consider a gas of Bose atoms with an internal en- 
ergy level structure as depicted in Fig. 1(a). The atoms 
populate the state \g) and form a Bose-Einstein conden- 
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sate characterized by the field-operator ip g (r) = p + 
Aipg (r), where po = Nq/Vq is the density of the (homo- 
geneous) condensate JlS| and the field operator A-tp g (r) 
describes small deviations from an ideal condensate due 
to non-condensate atoms. The Bose field-operators that 
account for the other two electronic states are expanded 
in terms of plane waves as 



^(r) = V /2 I>.ke 



(1) 



Polaritons are the normal modes of a combined system 
of radiation and matter and are a familiar subject in solid 



where 6kj annihilates an atom with internal state \j) 
(j = 9: e ) an d wave-vector k. A classical control field, 



with wave- vector k c and angular frequency w c , couples 
states q) and |e). In the presence of this field, the Hamil- 
tonian for the excited atoms reads 

-ffatom = ft J^ [w e ,k&;, k & e,k + (^g,k-k c + UcW^bq,* 
k 

+ ^ & g,k-kA k + fi C & e,k 6 «,k-k c ] , (2) 

where cjj.k = Wj + Tik 2 /2M (j — q, e) are the atomic 
excitation energies (in units of h), and fi c is the Rabi 
frequency determining the magnitude of the excited-state 
splitting. To remove an explicit time dependence due 
to the control field we have adopted a rotating frame 
(frequency lu c ) for the state \q). 

Adding now the interaction with a quantized probe 
field, the Hamiltonian for the combined system of radia- 
tion and matter is written 



H = H a 



H T 



;ad 



-^rad — atom T" ^o 



Here i? ra d is the radiative Hamiltonian, 



-ffrad = 2J ftc^(a k Ok + 1/2), 



(3) 



(4) 



p , in the probe pulse wave-packet; the operator 
atom = — ^o l Jd(r)p(r)d 3 r describes the interac- 



where a^ is the annihilation operator for a probe pho- 
ton (polarized along the dipole moment \x of the g — > e 
transition), and the summation is over all wave- vectors, 
k = k 

tion between the probe field and the atoms, where d(r) 
and p(r) = pijj e (r) ip g (r) + h.c. are the electric displace- 
ment and polarization field operators, respectively; the 
last term, fl" C ont = (x/2eo)J p(r)p(r)d 3 r, represents the 
contact interaction which appears in the multipolar for- 
mulation of QED |l9| ^2J , commonly adopted with x = 1 . 
We choose x = 2/3 in the contact interaction |E3|. This 
yields the correct local field corrections in the refractive 
index [Eq. (§]. 

As a first step we take into account only coherent op- 
tical processes in which the absorption of a photon pro- 
motes a condensate atom to the excited electronic state 
from which it returns, via photon emission, to the con- 
densate. For this purpose, we make the replacement 
tp g (r) — » p Q ' in the polarization field p(r). In doing 
this, incoherent processes, such as spontaneous emission 
to modes other than that of the absorbed photon, are 
disregarded. Spontaneous emission will be analyzed in a 
second step later on. With the replacement, one has 
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(5) 



and 



Hcont = XJ2 f/(^,-k + ^ k )(& e ,k + ^_ k ), (6) 



where the transition dipole moment is chosen to be real. 
The complete Hamiltonian, with the simplified terms 
(||) and (||), has been diagonalized by a Bogoliubov-type 
transformation, which mixes the radiation and matter 
modes to yield the polaritons. We thus obtain, adapting 
the methods of Juzeliunas [01 , 



h = n J2 J2 < p U p ™* + const 



(7) 



k m— 1 



where m = 1, 2, 3 labels the polariton dispersion branches 
(Fig. 2) and the eigen-frequencies w^ — w are to be de- 
termined from the equation to = ck/n, where n = n (u, k) 
is the refractive index, with 



n 2 = 



1 + xap/eo 



(x = 2/3) : 



1 - (1 - x) ap/e 
a is the atomic polarizability, given for |Aa->| <w by 
[i 2 Aw 



h Aw 2 + /3Auj - n 2 '' 



(8) 



(9) 



where Aw = lu — w c — w^k-ke is the detuning from the 
two-photon resonance and (3 = w g k-k c + u c — uj e ^ is the 
control laser frequency mismatch. 
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FIG. 2. Polariton dispersion branches (solid lines) in an 
EIT medium with /3 — 0. The dashed lines show the un- 
coupled modes in the limit of vanishing density. The dotted 
curve is calculated from the semiclassical susceptibility and 
includes the excited state linewidth. The parameters in (b) 
are taken from the experiment of Liu et al. [6] ; note the ex- 
panded frequency scale around the EIT resonance where the 
free-space photon branch is almost a vertical line. In (a) the 
density is reduced to give v g /c = 10 -1 so that the curves can 
be plotted using a single frequency scale. 



The polariton modes are defined by the Bose operators 
-Pm.k- These are determined explicitly to be (|Aw| <C to) 

PmM = P lnM + Mrifck-k. + »« (Aw/fi c ) &e,k, (10) 

where u = u m ,k = (1 — nv g/ c ) [1+ (Aw/Sl,-) 2 ] -1 / 2 and 



D (rad) 
m,k 



= K/c) 1/2 [(n + l)a k + ( ? i - l)aL k ]/2 (11) 



is the radiative component; here v fl = c (n + ojdn/duo) 
is the radiative group velocity, and both n = rv- m ' and 
w g = Vg are to be calculated at u> = w k . 

In contrast to the polaritons considered in [|15| [L7| , 
the operator P2,k represents an eigen-mode of the to- 
tal Hamiltonian, not only at EIT resonance (Aw = 0), 
but for the entire branch of slow polaritons (m = 2). 
The systematic k-space formulation we have presented 
facilitates a discussion of the frequency and wave-vector 
dependence of the radiation and matter fields. 

The polariton dispersion branches are plotted in the 
vicinity of the EIT resonance in Fig. 2. Figure 2(b), in 
particular, adopts the parameters of the recent experi- 
ment, in which v g /c « 2fl" 2 heo/cJolJ 2 p ~ 10~ 7 M. To aid 
comparisons with the EIT literature, where the refrac- 
tive index is usually plotted, we also plot the dispersion 
curve (dotted line) given by the semiclassical susceptibil- 
ity. The susceptibility includes the effect of the excited 
state lifetime and hence the dispersion curve continuously 
connects the three polariton branches. 

In our theory the polaritons acquire a finite lifetime 
through the excited-state contribution to Eq. (|10|). The 
decay rate is T = u 2 (Aw/J7 c ) To, where Tq is the rate of 
atomic radiative decay and u rs 1 provided |Aw| <C Q c - 
Away from the EIT resonance, the polariton decays spon- 
taneously at this rate into a photon and a translationally 
excited ground-state atom. In the experiment [pi the 
probe pulse excites a frequency band |Aw|/J7 c ~ 0.02. 
There is thus significant loss (~ 25%) due to polariton de- 
cay during the measured 11.8 /is propagation time in the 
medium. Some broadening of the probe pulse is expected 
from this as the decay rate is frequency dependent. 

The contact interaction has no discernable effect on 
the slow polariton branch of Fig. 2. At the density cur- 
rently used H its influence is at the level of less than 
1%. With a ten-fold increase in density, however, lo- 
cal field effects would no longer be negligible. The dis- 
tinction between the radiative and full group velocity 
du/dk = d(Au)/dk + duj q ^ k Jdk ~ v g + h(k - k c )/M 
is also unimportant in current experiments. The atomic 
velocity, h(k — k c ) /M, is much smaller than v g even for 
counter-propagating control and probe lasers (typically 
it is of the order of a few cm/s). 

We consider now the time evolution of a wave-packet of 
EIT polaritons (a probe pulse in the medium) influenced 
by the following sudden perturbations: (i) at t = t\ the 
original control beam is switched off, (ii) after some delay, 



at t = ti > 1 1, the control beam is again turned on with a 
new Rabi frequency fl' c and wavevector k^.. In contrast to 
what has been previously considered (|,|7],[l6|,|l7j , k.' c does 
not necessarily coincide with the original wavevector k c . 
If kj. = — k c , for instance, the regenerated probe pulse 
can move backwards, as we will show. 

At t < t\ the wave-packet is described by the state- 
vector 



I*) = |{«k,t}) 



Y\ exp(a k:t P 2 f k - h.c.) |vac) 



(12) 



where |vac) is the vacuum state- vector (P2,k |vac) = 0), 
|{c*k,t}) is a many-mode coherent state, and the ampli- 

tude ak,t = ak.ti exp[— zw k (t — t{)\. At t = t\ the con- 
trol laser is suddenly switched off. Equation (O) then 
represents the initial condition for the subsequent time 
evolution, giving for t\ <t <t%, 



|i)=JJexp| 



lUQfk 



tbl 



q.k— k c 



h.c.) exp(Sk)|vac). (13) 



The first exponent describes the magnetic excitations, 
which are now decoupled from the photons so that (up to 
a phase) a^.t = oik.ti exp [—ioj qt k-k e (t — ti)]. The opera- 
tor exp(Sk) accounts for the small radiative (~3x 10~ 4 ) 
and excited state (~ 0.02) contributions to the slow po- 
lariton amplitude. These excitations are subsequently 
converted into spontaneously emitted photons and may 
be omitted for t — t\ larger than the radiative lifetime. 
It should be emphasized that while the radiative contri- 
bution plays an essential role in giving the polariton its 
velocity while the control laser is "on," it accounts for 
a very small fraction of the quasiparticle number (when 
Vg/c <C 1) and may simply be discarded to stop the po- 
lariton. 

At t = ti the new control beam is applied. Expanding 
the magnetic operator ib' kk = ib k ,_ k / in terms of the 

new polariton operators P^ k , (m = 1,2,3), the state- 
vector evolves for t > ti as 

\t) = JJexp(Wa k / )t Pa| v - h.c.) exp(S" k ,)|vac), (14) 



with a;,, ,. = a kjt2 exp[— iuj^ 



*k'* 



/( , 2) (i - t 2 )\. The operator S v 



accounts for the other polariton modes (m =1,3) which 
are subsequently converted into spontaneously emitted 
photons. Here k' = k — k c + k^, is the wave- vector of the 
regenerated polariton and the quantities u', -P^k'; an( ^ 



Am) 



are defined in the obvious way 



Such behavior may be viewed as a kind of time-delayed 
four- wave mixing (photon echo |25| ) involving the orig- 
inal probe pulse, k, two control beams, k c and k^,, and 
the regenerated probe pulse, k', the wave- vectors satis- 
fying the phase matching condition k — k c = k' — \s.' c . 
The distinctive feature is the involvement of slow EIT 
polaritons. Consequently very high conversion efheency 



can be reached (uv! rs 1) if both the original pulse and 
the regenerated polaritons are in the vicinity of the EIT 
resonance. Note that there is no need for a smooth (adi- 
abatic) "turn-off" and "turn-on" of the control laser to 
map the slow polariton into and out of a magnetic ex- 
citation. Even if the switching on and off are instanta- 
neous, the reconstruction is almost perfect when uu' ~ 1. 
The actual storage and retrieval of the probe pulse takes 
place at the medium boundaries where vacuum photons 
are converted into slow polaritons and vica versa. 

We consider, finally, some specific situations. If the 
wavevector of the control laser does not change (k' c = k c ), 
one arrives at degenerate four wave mixing in which the 
regenerated probe photon has the same wavevector as 
the original (k' = k). This is the situation investigated 
in recent work P,[^, [16|jl7| l . On the other hand, with co- 
propagating control and probe beams one can change the 
direction of the probe pulse by changing the direction 
of the control beam |ej. If k' c = -k c (27), then k' = 
— k — 2 (k c — k), and the reverted polariton experiences 
a frequency shift Su>' « 2v'(k c — k). Under the conditions 
of the recent experiment || one obtains 6u>' rs 0.26 kHz x 
2-7T, which is very small compared to the coupling Rabi 
frequency (and also the spectral width of the probe pulse) 
and the reverted polariton remains in the EIT region. For 
counter-propagating beams, the wave-vector of the new 
polariton is k' = k + 2k' c w 3k; the polariton becomes 
situated close to the upper band edge where there is a 
rate of radiative decay on the order of Tq (unless il' c is 
extremely small) [E8| . 

We have formulated a theory of slow polaritons in 
atomic gases and applied it to the slowing down, stor- 
ing, and redirecting of laser pulses in an EIT medium. 
The polariton modes have been determined through a full 
diagonalization of the dissipationless Hamiltonian, and 
lifetimes introduced as a secondary step. With detuning 
included various four-wave mixing possibilities were an- 
alyzed. The possibility of reverting a stopped polariton 
was demonstrated. 
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